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Abstract 

It is known that the Entropy Power Inequahty (EPI) always holds if the 
random variables have density. Not much work has been done to identify dis- 
crete distributions for which the inequality holds with the differential entropy 
replaced by the discrete entropy. Harremoes and Vignat showed that it holds 
for the pair {B{m,p), B{n,p)), m, n S N, (where B{n,p) is a Binomial distri- 
bution with n trials each with success probability p) for p = 0.5 . In this paper, 
we considerably expand the set of Binomial distributions for which the inequal- 
ity holds and, in particular, identify no{p) such that for all m, n > nQ{p), the 
EPI holds for {B{m,p), B{n,p)). We further show that the EPI holds for the 
discrete random variables that can be expressed as the sum of n independent 
identical distributed (IID) discrete random variables for large n. 

1 Introduction 

The Entropy Power Inequality 

^2h(X+Y) y ^2h(X) _^ g2h(Y) ^-^-^ 

holds for independent random variables X and Y with densities, where h{-) is the 
differential entropy. It was first stated by Shannon in Ref. and the proof was 
given by Stam and Blachman 0. See also Refs. i, i, 0, 0, i, B- 

This inequality is, in general, not true for discrete distributions where the dif- 
ferential entropy is replaced by the discrete entropy. For some special cases (binary 
random variables with modulo 2 addition), results have been provided by Shamai and 



Wyner in Ref. [10 
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More recently, Harremoes and Vignat have shown that this inequahty will hold if 



X and Y are B{n, 1/2) and B{m, 1/2) respectively for all m,n 11|. Significantly, the 



convolution operation to get the distribution of X + F is performed over the usual 
addition over reals and not over finite fields. 

Recently, another approach has been expounded by Harremoes et. al. 12| and by 
Johnson and Yu 13|], wherein they interpret Renyi's thinning operation on a discrete 
random variable as a discrete analog of the scaling operation for continuous random 
variables. They provide inequalities for the convolutions of thinned discrete random 
variables that can be interpreted as the discrete analogs of the ones for the continuous 
case. 

In this paper, we take a re-look at the Harremoes and Vignat HI] result for the 
Binomial family and extend it for all p G (0, 1). We show that there always exists an 
no{p) that is a function of p, such that for all m,n> no{p), 



2H[B{m+n,p)] y ^2H[B{m,p)] _^ ^2H[B{n,p)] 



+ e^ 



(2) 



where H{-) is the discrete entropy. The result in Ref. [ll[ is a special case of our result 
since we obtain no(0.5) = 7 and it can be checked numerically by using a sufficient 
condition that the inequality holds for 1 < m, n < 6. 

We then extend our results for the family of discrete random variables that can be 
written as the sum of n IID random variables and show that for large n, EPI holds. 

We also look at the semi-asymptotic case for the distributions B{m^p) with m 
small and B{n,p) with n large. We show that even when n is large, there may exist 
some m such that EPI may not hold. 

Lastly, we show that how the EPI for the discrete case can be interpreted as an 
improvement to the bounds given by Tulino and Verdu for special cases [3]. 



2 EPI for the Binomial distribution 

Our aim is to have an estimate on the threshold no{p) such that 

^2H[B{'m+n,p)] y ^2H[B{m,p)] ^2H[B(n,p)] ^ j^g^j 

holds for all 771,11 > no(p). 

It is observed that no{p) depends on the skewness of the associated Bernoulli 
distribution. Skewness of a probability distribution is defined as K3/a/k| where K2 
and ^3 are respectively the second and third cumulants of the Bernoulli distribution 
-8(1, p), and it turns out to be {2p — 1)/ \/p{l — p)- Let 

M.) ^ (4) 

P{1 -p) 

We find an expression for rtQip) that depends on oj{p). The following theorem, known 
as Taylor's theorem, will be useful for this purpose (see for example p. 110 in Ref. 

3)- 
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Theorem 1 (Taylor). Suppose f is a real function on [a,b], n G N, the {n — l)th 
derivative of f denoted by /("~^) is continuous on [a, 6], and f^'^\t) exists for all t G 
(a, 6). Let a, f3 be distinct points of [a,b], then there exists a point y between a and 
P such that 



n-l 



= /(«) + E ^^^^ - + ^^(/^ - ")'^- (5) 

k=l 

For < p < 1, let H{p) denote the discrete entropy of a BernouUi distribution 
with probabihty of success p, that is, H{p) = — plog(p) — {1 — p) log(l — p). We shaU 
use the natural logarithm throughout this paper. Note that we earlier defined H{-) 
to be the discrete entropy of a discrete random variable. The definition to be used 
would be amply clear from the context in what follows. Let 

H{x) ^ H{p) - H{x), a;G(0,l). (6) 

Note that H{x) satisfies the assumptions in the Theorem [1] in x G (0, 1). Therefore, 
we can write 



H{x) = H{p) + J2 ^^(^ - p) + ^M^(^ - P)^ 

k=l ' ^' 

for some Xi G {x,p). Note that H{p) = and 

riog(x) - log(l - x), ifA; = l, 



(7) 



k\ 



1 



(l-x)~^'^^'^ + (-l)V(^-i) 



fc(fc-i) 

For even k, F^''\x) > for all x G (0, 1), and hence, 

21+1 



iik>2. 



H{x)>Y^F'^^\p){x-pf (9) 



fc=i 



for all X G (0, 1) and any non-negative integer /. The following useful identity would 
be employed at times 

iog(2)-i/(rt=i:^-j^(p--) , (10) 

Let P = {pi} and Q = {qi} be two probability measures over a finite alphabet A. 
Let C(p)(P,g) and aIJ'\p,Q) be measures of discrimination defined as 

C^^\P,Q) ^ pD{P II M) + qD{Q II M), (11) 



where 



A 



M ^ pP + qQ, 



1 — p. 



These quantities are generahzed capacitory discrimination and triangular discrimina- 
tion of order v respectively that were introduced by Tops0e [Isf . 

The following theorem relates C^^\P, Q) with aI^\p, Q) and would be used later 



to derive an expression for no{p). It generalizes Theorem 1 in Ref. 15 



Theorem 2. Let P and Q be two distributions over the alphabet A and < p < 1. 
Then 



!/ = l ^ ' 



Proof. Let 



We have 



rrii = ppi + qqi, = \2ppi - | , ki 



mi 



J_ _ \pPi - qqi 
ki pPi + qqi 



and 0<l/ki<l. We have 



ieA 



C^'\P,Q) = pJ2P^logi^)+qJ2Q^^ogi^ 



Pi log 



ieA 



PPi 



i&A 

q^qi loj 

ieA 



qi_ 

m, 



a \ ^ mi + ej 

= 2^—^ log 



i£A 



E 

ieA 



qqi 

mi 
mi - e,; 



H{p) 



log 



mi 



- [log(2) - H{p)] 

+ h) log (l + + - 1) log (l - i 

A \ i / \ i 



[log(2) - H{p)] 



ieA 

ie^ i/=i 

°° Ai^)(p,g) 



log(2) - if ( 1 + -L 



[log(2) - H 



2 '2ik^ 

- [log(2) - 



2iy{2iy-l)k 

[log(2)-/7(p)], 



(13) 

(14) 
(15) 

(16) 
(17) 



(18) 

(19) 
(20) 

(21) 

(22) 



where a follows by taking two cases 2ppi > nii and 2ppi < nii, and b follows from 
(HOD. □ 

Let X^"^^ be a discrete random variable that can be written as 

X^*^) = Z, + Z2 + --- + Zn, (23) 

where Zj's are IID random variables. We note that when X^^'' is defined as above, 
we have X^"^ + = Let K = emix^"^)], We first use a lemma due to 



Harremoes and Vignat [11 



Lemma 1 (Harremoes and Vignat ll|). If Yn/n is increasing, then Yn is super- 
additive, i.e., Ym+n ^Ym + Yn. 

It is not difficult to show that this is a sufficient condition for the EPl to hold 



ll[. By the above lemma, the inequahty 



i7(X(-+i)) - i/(X(")) > - log f ) (24) 



2 V n 

is sufficient for EPl to hold. 
Let = B{n,p). We have 

Px("+i) {k + l)= pPxi^) (k) + qPxM (k + l). (25) 

Define a random variable X^"^ + 1 as 

{k + l) = Pxi,4k). (26) 

for all A; G {0, 1, ■ ■ ■ , n}. Hence, using i/(X(") + 1) = iy(X(")), we have 

+ (27) 
i/(X("+i)) = pi/(X(") + 1) + gi/(X(")) +pD(P^(„,+i II P;,(„+i)) + gD(P;,w II P^(„+i)) (28) 
= i/(X(")) +pD(P^(„)+i II Px(n+i)) + g/^(Pxw II Px("+i))- (29) 

Therefore, 

i/(X("+^)) = i7(X(")) + C(P)(P^(„)+i,Px("))- (30) 
We now derive the lower bound for C^^'^ (^Px(n)_^_i, Pxi-n,) ) . 

Lemma 2. For I E N, 

n+l 



C(P)(P^(„)+i,P^,„,) = J2H(^)Pxi'^M^), (31) 

C(^)(P;,(„),P;,(„)+l) > 5^P(''Hp)(r^+l)"Vl"^'\ (32) 



fc=l 
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where /i^"'' is the k-th central moment of B{n,p), i.e., 

n 



i=0 



Proof. Let P = X(") + 1 and Q = We have 



pp^ + qqi + 

22 - n - 1 

72 + 1 ' 



(33) 



(34) 

(35) 
(36) 



i=o ^ ^ 



2 \ s-^ I . n + 1 



n 



1=0 



(37) 



Pxin+r>{l). (38) 



Using Theorem [21 we have 



-[\og{2)-Hip)] 

n+l oo 

EE 



Px(.n + l){i) 



-)2u 



= E 



2i^(2z/- 1) + 1 2 
-[log(2)-i/(p)] 



2;^ 



(39) 



(40) 



n+l r 



i=0 



log(2) - 



P^(„+,)(2) + i/(p)-log(2) (41) 



n+l 



i=0 



n + l 



Px(n+i){i) 



n+l 

E^ 



i=0 



n + l 



Px(n+l){i), 



(42) 
(43) 



where 'a' follows by using (fTOj) and '6' follows by using (Q. To prove the lower bound, 
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we have 

n+l 2/+1 / . s fc 



i=0 fc=l ^ ^ 

2i+l n+l . . . fc 

fc=l j=0 ^ ^ 

k=l 

where 'a' holds for all nonnegative integers / using (Q. □ 

The following lemma shows that unlike the continuous case, EPI may not always 
hold. 

Lemma 3. For p ^ 0.5, EPI does not hold for all n. 
Proof. It suffices to show that 

^2HlBi2,p)] ^ ^2HlBil,p)] _^ g2H[B{l,p)] y (^^j^ 

with equality if and only if p = 0.5. In other words, we need to show that 

H[Bi2,p)]-H[B{l,p)]-Uog{2) <0 W p 0.5. (48) 

Using Lemma |2] and (ITOl) , we have 

H[Bi2,p)]-H[B{l,p)] = H{p)-2p{l-p)\ogi2), (49) 

H{p) < log(2) - 2 (p - 0.5)^ . (50) 

Therefore, 

H[B{2,p)]-H[B{l,p)]-^^ < i^-2(p-0.5)^-2p(l-p)log(2) (51) 

= 2 (p - 0.5)' [log(2) - 1] (52) 
< Oifp^O.5. (53) 

In other words, EPI holds for Binomial distributions B{n,p) for all n only if p = 
0.5. □ 

For the case m = 1 and n = 2, Fig. [1] shows the plot of 

/(m, n, p) ^ e2^[^('"+"'*')l - {e^^t^^'"'^)] + e'^I^^"'^)! } (54) 

as a function of p. Note that EPI is satisfied for p close to 0.5, while EPI does not 
hold if p is close to or 1. 

This leads us to the question that for a given p, what should m, n be such that 
the EPI would hold. The main theorem of this section answers this question. 
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Figure 1: Plot of e^^^I^^^.p)] _ |e2//[B(2,p)] ^ g2/f[B(i,p)]| 

versus p. 

Theorem 3. 

+ > -log f j V n>no(p). 



(55) 



Several candidates of nQ{p) are possible such as no{p) = 4.44 uj{p) + 7 and no{p) 
u{pf + 2.34 uj{p) + 7. 



□ 



Proof. See Appendix Rl 
2.1 Lower bound for the entropy of the Binomial distribution 



Unlike the asymptotic expansion of H[B{n,p)] given in Ref. 16|], we give non- 
asymptotic lower bound to it. Let 



21+1 



k=l 



We have 



where X*^"^ = B(n,p). 



i7(x(^))-i7(x(^-^))>rKj), 



(56) 



(57) 



Using Faa di Bruno's formula 17|, we have 

Z^if^ = y ■ u n- ^' ■ u n- < ■ ■ ■ < , (58) 

where Kg is the g-th cumulant of the BernouUi distribution and 

hgi + i2g2 H h isgs = k. (59) 

The summation is over all such partitions of k. We have 

21+1 



fe=i 

21 

= Y.r'^ciw), (61) 



w=l 



where 



c{w) = y ———————nl • ■ ■ < F^^^ (p) (62) 

where the summation is over all such (ii, 12, - ■ ■ , is] gii g2, - ■ ■ , gs) such that 

ii{gi - 1) + i2{g2 - 1) H \- is{gs - I) =w. (63) 

Now 

n 

i7(X(")) = H{X^'^^) + J2[HiX^'^)-H{X^'~^^)] (64) 

n 

> ^r,(j), since i/(X(o)) = (65) 

n 21 

= EE^'"'"^^ (66) 

j = l UI = 1 

= E^(^)E^'"'"- (67) 

«)=1 J = l 

Note that y^""., are the Generalized Harmonic Numbers (see for example Ref. 

0). 

As an example, we compute the lower bound for if (X*^"^) for I = 1. We have 
c(l) = K2F^'^\p) and c(2) = ?)k\F^^\p) + fi:3F'^^^(p). The first and second cumulants 
of Bernoulli distribution is given by K2 = p{l — p) and K3 = p{l — p){l — 2p). This 
gives c(l) = 1/2 and c(2) = [1 — J9(l — j9)]/[12j9(l —p)] and we get the lower bound as 
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3 EPI for the sum of IID 

We showed in the previous section that EPI holds for the pair {B{n,p), B{m,p)) for 
all m,n > tiqIp). The question naturally arises whether EPI holds for all such discrete 
random variables that can be expressed as sum of IID random variables. Let X^'^^ be 
a discrete random variable such that 

X(") ^Xi+X2 + --- + Xn, (69) 

where Xj's are IID random variables and cx^ is the variance of Xi. We shall use the 
asymptotic expansion due to Knessl [16 . 



Lemma 4 (Knessl [16|). For a random variable X^'^\ as defined above, having finite 
moments, we have as n oo, 

,,(„) 4 „(X-^) - i |og(2™n.') -^i + 53 JL, (70) 

1=1 

where Kj is the jth cumulant of of Xi. If = K4 = ■ ■ ■ = Kn = but Kn+i 7^ 0, then 

^ ' l=N-l 

Note that the leading term in the asymptotic expansion is always negative. We 
also note using Lemma H] that as n — 00, 

i/(X(")) < llog(27rena2). (72) 

To see this, we invoke the definition of the asymptotic series to get 

/ \ 1-N , Pn-1 , f ^ \ /^7Q^ 

Prom the definition of the "little-oh" notation, we know that given any e > 0, there 
exists a L(e) > such that for all n > L{e), 

9W 2(iV+l)!a27V+2^ + ■ 

Choosing n large enough, we get the desired result. 
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3.1 Asymptotic case 

We first consider the case of tlie pair wlien botli m, n are large and liave 

the following result. 



Theorem 4. There exists a no € N such that 



(75) 



for all m,n > uq. 



Proof. We shall prove the sufficient condition for the EPI to hold (as per Lemma [T]) 
and show that 



H{X 



> - log 



1 



n + 1 



n 



(76) 



for n > Uq for some Uq G M. 

Let us take the first three terms in the above asymptotic series as 



Co 



n' 



k2 



where < ki < k2 < and Ci is some non-zero positive constant, and hence, 



9{n) + 

1 f^l 



1 



and given any e > 0, there exists a L(e) > such that for all n > L{e) 



9W + -^fc, 



Co C; 



3_ 



< e 



Therefore, we have the following inequality 



C2 _|_ C*3 



n 



k2 



< g{n) < —J- + -J- 



(77) 



(7J 



(79) 



^0) 



From inequality ( IHOj) . we can say by using the lower and upper bounds respectively 
for g{n + 1) and g{n) that. 



g{n + 1) - g{n) > Ci 



{n + iy 



c, 



1 



1 



+ 



_(n + l)^'2 n''\ 
Cs-e Cs + e 



_(n + 1)^'3 



n 



(81) 



From the above expression, we can clearly see that the first term is strictly positive 
and is O (l/n'^^"'"^) . The second and third terms (their signs are irrelevant) are of the 
order 0{l/n^^~^^) and 0{l/n^'-^) respectively. It is clear that there exists some positive 
integer no such that for all n > uq, first (positive) term will dominate and the other 
two terms will be negligible compared to the first and hence g{n+l) — g{n) > 0. □ 
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3.2 Semi-asymptotic case 

We now consider the pair where n — )■ oo and m is fixed. We already 

know from the previous result that the EPI holds when both m and n are large. 
We start by writing an asymptotic expansion of 



/(m, n) = e 



A 2H(X('"+")) _ 2_H'(X(")) _ 2_H'(X('")) 



^2) 



by using Knessl's result in LemmalHin which an asymptotic expansion for the entropy 
of X^"^ is derived as 



^ oo 

~ - log(27rena2) + ^ 



n 



Let 



We have 



2^ ^i+i ■ 



1=0 



(m+n) 



[27re(m + n)cr2]e5('"+"), 



(83) 



(85) 
(86) 



and we can rewrite 

/(m, n) = Wnema^ - e2^^(^""')l + 27re(m + n)a^ [e^(™+") _ l] 

- 27iena^ [e^(") - l] . 



^7) 



The first term in the above equation is a constant since it depends only on m and 
the second term can be written as 



^ f +o( 1 "ll 

27re(m + ?7,)cr <eL™+" ('"+")^ V(m+n)^;j _ i 



which can be expanded into 
27re(m + n)a'^ 



(3o 



+ 



2/3l + P'n 



m + n 2(m + ?7,)2 
Similarly, the third term can be written as 



+ 



[m + nY" 



-2'Kena 



+ o( - 



n 2?T,^ \ n 

Using the above two expressions, o [l/(m + n)^] = o (1/n^), and (E 

'2/3i + /3o'Y^ 1 

m + n n 



f{m,n) = [27rema2 - e2^^(^*'"')] + 2W 



we get 
n / 



(89) 



(90) 



(91) 
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where the terms 

1 

n 



o 



and 2TTea^ 



2(3, + (3i\f 1 



m + n n 



can be made arbitrarily small as n — oo. Therefore, for large enough n, we can 
see that the first term dominates over other terms and moreover, f{n,m) > if 
2Trema^ - e2^^(^''"') > 0. Therefore, 

if n — OO and 

< ilogpTremd^]. (93) 

It follows from f l72|) that the above inequality holds for sufficiently large m. For the 
Binomial distribution B{m,p), EPI will hold for all such p that satisfy 

H[B{m,p)] < ^log[27remp(l -p)]. (94) 
The above relation is not true for all p and m. 

4 Discussion and Conclusions 

We show that how our results can be used to improve a bound by Tulino and Verdu 
0] under special cases. 

4.1 Improvement on a bound by Tulino and Verdu 

Let Xj, i = l,2,...,n be discrete IID random variables and Zj, i = 1,2, ...,n, be IID 
random variables as well with Zi ~ A/'(0, cx^). Let 

n 

= Y,i^^ + Z^) (95) 

1=1 
n 

i=l 

Let 

D{Y) = 0.5 log(27re4) - /i(F), (97) 

where F is a random variable with density and variance Uy. Tulino and Verdii [3l 
interpreted D{Y) as the non-Gaussianess of the random variable X and showed that 
the non-Gaussianess increases by having more random variables, i.e., 

D(^(")) < Z^(^("-^)). (98) 
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Expanding it using ( 1971) . we get 

/.(^W) - h{S^--'^) > \ log (^^) . (99) 
We show that for sufficiently large n, this bound can be made tighter for small a, 



I.e., 



lim [/i(^(")) - /i(5("-i))] > log (^;^) • 



(100) 



Let 

/„ = /(X("); ^(")) = /i(^(")) - 0.5 log(27rencr'). (101) 
Note that using Lemma 6 in Ref. [l9| 

i7(xW) = lim(/„). (102) 

(J— >o 

We know that for sufficiently large n, the EPI holds and 



i7(X("))-i7(X(""i))>ilog(^^ 



(103) 



Therefore, 

Ito (4-4^.)>ilog(^) (104) 



for sufficiently large n. The above limit is due to [19[. On the other hand, 

/„ - = h{S^^^) - /i(^("-i)) - \ log (-^\ ■ (105) 

2 \n — I J 

Hence, 



lim [/i(5(")) - /i(5("-^))] > log f (106) 
o-^o \r2 — 1/ 

for sufficiently large n. Comparing (11061) with fl99|l . we note that our bound is tighter 
by a factor of 2. 

4.2 Conclusions 

In conclusion, we have expanded the set of pairs of Binomial distributions for which 
the EPI holds. We identified a threshold that is a function of the probability of 
success beyond which the EPI holds. We further show that EPI would hold for 
discrete random variables that can be written as sum of IID random variables. 

It would be interesting to know if C*'^^^(Pjf(„)_,_]^, Pj^C")) ^^'^^ = B{n,p) is a 
concave function in p. It would also be of interest to know that for a given p G (0, 0.5) 
if H[B{n + l,p)] — H[B{n,p)] — 0.51og(l + 1/n) would have a single zero crossing as 
a function of n when n increases from 1 to oo. 
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A Proof of Theorem [3 

We prove that 



H[B{n + l,p)] - H\B{n,p)] >-\og 



n + l 



n 



W n > rioij)). 



Using ( H6|l . we have 

H[B{n + l,p)]-H[B{n,p)] = C(p)(P^(„), P^^+i) 



k=l 



Let 



r = p — 1/2. 

We have the first seven central moments of B[n,p) as 



10n(3-64r2 + 208r^)], 



(n) 



32 



8(17- 240r2 + 720r% 



(107) 

(108) 
(109) 

(110) 



in) 

= 


^n(l -4r2), 


(111) 


(n) 

= 


— -nr(l — 4r^), 


(112) 


(n) 

= 


— n(l - 4r2)[-2 + 24r2 + 3n(l - Ar^)], 
16 


(113) 


(n) 

= 


~nr{l - 4r^)[-4 + 24r^ + 5n(l - 4r^)], 


(114) 


(n) 
f^6 = 


—72(1 - 4r2) [15^2(1 - Ar^ f + 16(1 - 30r2 + 120r^) 
32 





(115) 



nr{l - 4r^)[105n^(l - Ary - 14n(17 - 200r^ + 528r*), (116) 



(117) 



Let t = Lo{p) = 16rV(l -4r2) and hence, = t/[4(t + 4)]. Note that e [0, 1/4) 
and t G [0, oo). The above seven central moments contain only even powers of r and 
hence, can be written as a function of t. 

We upper bound the right hand side of (11071) as 



log 



n + l 



n 



^ 1 1 1 



(118) 



Define 



k=l 



t 1 

+ 



4(t + 4) 2 



1 1 1 

n 271^ 3?7,2 



(119) 
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Proving f llOTp is equivalent to showing that /(n, t) > V n > rioip). Simplifying 



f{n,t) 




(-2989t - 721t^ - 3339t^ - 315)n^ + 
(721t - 546t^ - 826 + 371^^ - 1568t^)n^ + 
(-135^2 - 157t^ - - 826 - 90t - 66t^)n' 

630^2 - 315n - 70 . 



3 



(120) 



Define 



^(n,t) = 420(n + l)V/(n,t). 



(121) 



A simple but elaborate calculation yields (7(4. 44t + 7 + m,t) ^ 35tm'' + (1122. 8t^ + 
2030t + 70)m'' + (14700.90t3 + 52210. 20^^ + 48120.80t + 2625)m5 + (l.Olt^ + 5.32^^ + 
9.57t2+6.06t+0.40)105m^ + (3.85t5 + 26.94t^ + 72.32t3+88.61t2+43.61t+3.02)10^m3 + 
(7.76^6 + 68.23t5 + 247.042t^ + 456.97t3 + 433.17^^ + 176.77t + 11.80)10 W + (6.47t^ + 
70.91t6 + 338.88t^ + 880.98t^ + 1297.85t=^ + 1030.51t2 + 361.59t + 20.14)10^m+ (0.15^8 + 
56.29t^+709.80t6+3485.03t^+8728.40t^+11955.74t3+8613.06t2+2628.77t+64.15)10^ 

Note that all the coefficients are positive and hence, /(4.44t + 7 + m, t) > for 
all m > or f{n, t) > for all n > 4.44t + 7. A more careful choice would yield all 
coefficients to be positive for n > 4.438t + 7. Yet another choice that would yield 
all coefficients as positive would be n > + 2.34t + 7. Note that this choice would 
be better for < t < 2.1 and, in particular, for t = 1, the first choice yields (after 
constraining n to be a natural number) n > 12 while the second one yields n > 11. 

Further refinements are also possible. For example, the expansion of /(7 + 
m, t/[4(l + t)]) yields positive coefficients again. Such a choice constrains uj{p) G 
(0, 1/4) and we get no{p) = 7. 
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